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Abstract 



A complete set of QCD sum rules for the magnetic moments of decuplet 
baryons are derived using the external field method. They are analyzed thor- 
oughly using a Monte-Carlo based procedure. Valid sum rules are identified 
under the criteria of OPE convergence and ground state dominance and their 
predictions are obtained. The performances of these sum rules are further 
compared and a favorable sum rule is designated for each member. Correla- 
tions between the input and the output parameters are examined and large 
sensitivities to the quark condensate magnetic susceptibility \ are found. Us- 
ing realistic estimates of the QCD input parameters, the uncertainties on 
the magnetic moments are found relatively large and they can be attributed 
mostly to the poorly-known \- It is shown that the accuracy can be improved 
to the 30% level, provided the uncertainties in the QCD input parameters 
can be determined to the 10% level. The computed magnetic moments are 
consistent with existing data. Comparisons with other calculations are made. 
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I. INTRODUCTION 



The QCD sum rule method ]l[ has proven a powerful tool in revealing the deep con- 
nection between hadron phenomenology and QCD vacuum structure via a few condensate 
parameters. The method has been successfully applied to a variety of problems to gain a 
field-theoretical understanding into the structure of hadrons. Calculations of the nucleon 
magnetic moments in the approach were first carried out in Refs. and ||. They were 
later refined and extended to the entire baryon octet in Refs. 0-0]. On the other hand, 
the magnetic moments of decuplet baryons were less well studied within the same approach. 
There were previous, unpublished reports in Ref. [§] on A ++ and Q~ magnetic moments. 
The magnetic form factor of A ++ in the low Q 2 region was calculated based on a rather 
different technique 0. In recent years, the magnetic moment of Q~ has been measured 
with remarkable accuracy [I0|: fi n ~ = —2.02 ± 0.05 fi N . The magnetic moment of A ++ 
has also been extracted from pion bremsstrahlung [11]: fi A ++ = 4.5 ± 1.0 fi N . In an earlier 
work []12| , the magnetic moment of A extracted from ir~p bremsstrahlung data was found 



to be consistent with /i A o = 0. The experimental information provides new incentives for 
theoretical scrutiny of these observables. 

In this work, we present a systematic, independent calculation of the magnetic moments 
for the entire decuplet family in the QCD sum rule approach. The goal is two-fold. First, we 
want to find out if the approach can be successfully applied to these observables by carrying 
out an explicit calculation. Second, we want to achieve some realistic understanding of the 
uncertainties involved in such a determination by employing a Monte-Carlo based analysis 
procedure. This would help us assess the limitations and find ways for improvements. 

We will show that both goals are achieved in this work. The entire calculation is more 
challenging than the octet case due to the more complex spin structure of spin-3/2 particles. 
One has to overcome enormous amount of algebra to arrive at the final results. But con- 
ceptually it presents no apparent difficulties. Particular attention is paid to the complete 
treatment of the phenomenological representation, which leads to the isolation of the tensor 
structures from which the QCD sum rules for the magnetic moments can be constructed. 
Flavor symmetry breakings in the strange quark are treated consistently across the decuplet 
family. The success also hinges upon a new analysis of the two-point functions ||13|| , which 
provides more accurately determined current couplings for normalization. Part of the results 
on A ++ and Q~ have been communicated in a letter Hl4 . 



Magnetic moments of decuplet baryons have also been studied in various other meth 

chiral perturbation theory 



ods, including lattice QCD [|T 

17] , non-relativistic quark model relativistic quark models |p"9|-|2 
chiral bag model [ 25J , cloudy bag model 



ism 



soliton model 24 



Bethe-Salpeter formal- 
chiral quark- 
A 



Skyrme model |27 

comparison will be made with some of the calculations and with existing data. 

Sec. H deals with the derivation of the QCD sum rules. Sec. [II J discusses the Monte 



Carlo analysis procedure. Sec. |TV| gives the results and discussions. Sec. [V] contains the 
conclusions. The Appendix collects the QCD sum rules derived. 
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II. METHOD 



Consider the time-ordered two-point correlation function in the QCD vacuum in the 
presence of a constant background electromagnetic field F^: 

IMp) = i J d 4 x e**<0 | T{ r] a (x) fj p (0) } | 0) F , (1) 

where r] a is the interpolating field for the propagating baryon. The subscript F means that 
the correlation function is to be evaluated with an electromagnetic interaction term added 
to the QCD Lagrangian: 

d = -A„J», (2) 

where is the external electromagnetic potential and J M = tqCpi^q the quark electromag- 
netic current. 

Since the external field can be made arbitrarily small, one can expand the correlation 
function 

n a ,(p) = nS(p) + nS(p) + ---. (3) 

Here n^ j (p) is the correlation function in the absence of the field, and gives rise to the mass 
sum rules of the baryons. The magnetic moments will be extracted from the QCD sum rules 
obtained from the linear response function H^(p). 

The action of the external electromagnetic field is two-fold: it couples directly to the 
quarks in the baryon interpolating fields, and it also polarizes the QCD vacuum. The latter 
can be described by introducing new parameters called vacuum susceptibilities. 

The interpolating field is constructed from quark fields, and has the quantum numbers 
of the baryon under consideration. We use the following interpolating fields for the baryon 
decuplet family: 
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(4) 
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Here implicit function forms r](x) and q(x) (q=u,d,s) are assumed. C is the charge con- 
jugation operator. The superscript T means transpose. The indices a, b and c are color 
indices running from one to three. The antisymmetric tensor e abc ensures the three quarks 
form a color singlet state. The normalization factors are chosen so that correlation func- 
tions of these interpolating fields coincide with each other under SU(3)-flavor symmetry (see 
Eqs. © to (H)). 

The interpolating field excites (or annihilates) the ground state as well as the excited 
states of the baryon from the QCD vacuum. The ability of a interpolating field to annihilate 
the ground state baryon into the QCD vacuum is described by a phenomenological parameter 
\b (called current coupling or pole residue), defined by the overlap 

(0\r] a \Bps) = X B u a (p, s), (5) 



where u a is the Rarita-Schwinger spin- vector [28 



A. Phenomenological Representation 

On the hadronic level, let us consider the linear response defined by 



IlS(p) = i jd A x e^(OMx) 



i d 4 y AMJ»(y) 



%(0)|0). (6) 



After inserting two complete sets of physical intermediate states, it becomes 

d A k d*k' 

4 

BB 



ng(P) -j d ' X j (27r) 4 (27r) 4 g ; £ k 2 _ M 2 _ le k ,2 _ ^ _ le 

e^( 2/ )(0|r /a (x)|A; S )(^|^( Z/ )|A;V)(fcV|r ?/3 (0)|0). (7) 



QCD sum rule calculations are most conveniently done in the fixed-point gauge. For elec- 
tromagnetic field, it is defined by x^A^i^x) = 0. In this gauge, the electromagnetic potential 
is given by 

AM = -i (8) 



The electromagnetic vertex of spin-3/2 baryons is defined by the current matrix element [15 



(ks\J^0)\k's')=u a (k,s) O a ^(P,q)up(k',s'). (9) 
The Lorentz covariant tensor 
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where P = k + k! and q = k — k', satisfies the standard requirements of invariance under 
time reversal, parity, G parity, and gauge transformations. The parameters ai, a 2 , c\ and c 2 
are independent covariant vertex functions. They are related to the multipole form factors 
by 

G E0 {q 2 ) = (1 + §t) [01 + (1 + r)a 2 ] - |r(l + r) [c x + (1 + r)c 2 ] 
G E 2{q 2 ) = [ai + (1 + r)a 2 ] - 1(1 + r) [c x + (1 + r)c 2 ] 

G M i(g 2 ) = (l + |r)a 1 -|r(l + r)c 1 ^ 
GW*? 2 ) = ai - 1(1 +r)ci. 



where r = —q 2 / (2Mb) 2 (> 0). They are referred to as charge (E0), electric quadrupole (E2), 
magnetic dipole (Ml), and magnetic octupole (M3) form factors. The magnetic moment 
is related to the magnetic dipole form factor Gmi(? 2 ) at zero momentum transfer. From 
Eq. (|TTD , it is clear that 

G M1 (0) = a 1 = fi B , (12) 



where the magnetic moment fi B is in units of particle's natural magneton: eft/(2cM#). So 
the goal is to isolate terms in Eq. ([?]) that involve only a\. 

The ground state contribution to Eq. (0) can be written as 
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(13) 



In arriving at Eq. (|T3|), we have used a number of steps: the translation invariance on rj a (x) 
and J^{y)i a change of variable from k! to q, the relation 

/ oV^V = -*(2tt) 4 ^(q), (14) 
integration by parts, and the Rarita-Swinger spin sum |28[1 

f \~ ( \ C I Jl/T \l 1 2 PaP(3 .Pal/3 ~Pf3la\ n 

2^u a (p, s)u^(p, s) = -(p + M B ) \ g a/ 3 - -7^ - + I , (15) 

with normalization u a u a = 2Mb- The hat notation denotes p = p a j a . 

Direct evaluation of Eq. ( |I~3"D leads to numerous tensor structures, not all of them are 
independent of each other. The dependencies can be removed by ordering the gamma ma- 
trices in a specific order. Here we choose to order in Plalixlulp- After a lengthy calculation, 
18 tensor structures which involve only a\ are isolated. They can be organized as 
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(16) 



The tensor structures associated with WEj have odd number of gamma matrices, while those 
associated with WOj have even number of gamma matrices. Apart from a common factor 
i\ 2 B H B /{p 2 — Ml) 2 , the invariant functions are given by 
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(17) 
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In addition to the ground state contribution, there exist also excited state contributions. 
For a generic invariant function, the pole structure has the form 
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(18) 



where Cb^b* are constants. The first term is the ground state double pole which contains 
the desired magnetic moment of the baryon, the second term represents the non-diagonal 
transitions between the ground state and the excited states caused by the external field, and 
the ellipses represent pure excited state contributions. Upon Borel transform, one has 
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Ml 



(1 



-{Ml, -M|)/M 2N 



+ 



(19) 



We see that the transitions give rise to a contribution that is not exponentially suppressed 
relative to the ground state. This is a general feature of the external-field technique. The 
strength of such transitions at each structure is a priori unknown and is an additional 
source of contamination in the determination of /j,b not found in mass sum rules. The 
usual treatment of the transitions is to approximate the quantity in the square brackets 
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by a constant, which is to be extracted from the sum rule along with the ground state 
property of interest. Inclusion of such contributions is necessary for the correct extraction of 
the magnetic moments. The pure excited state contributions are exponentially suppressed 
relative to the ground state and can be modeled in the usual way by introducing a continuum 
model and threshold parameter. 



B. Calculation of the QCD Side 



On the quark level, one evaluates the correlation function in Eq. ([!]) using Operator 
Product Expansion (OPE). The calculation is most readily done in coordinate space. To 
arrive at the final sum rules, one needs a subsequent Fourier transform, followed by a Borel 
transform. 

We decide to carry out four separate calculations for f2~(sss), S* + (uus), S*°(uss), and 
S*°(uds). They have distinct strange quark content, which requires special treatment. The 
QCD sum rules for other members can be obtained by appropriate substitutions in those 
for these four members. 



The master formula, which is obtained from contracting out the quark pairs in the 
correlation function, is given by, for Q~: 



(0\T{ t]^~ (a?) f)*p~ (0) } | 0) j 



2e abc e a'b'c' | S aa' Tl j p CSf C^Sf +2Sf^CSf Cj a Sf}, 



(20) 



for E* 



(Gin r^ro^n )}!^ 



ST Tr 



^CS b u b c la st 



T , 



+2S? i P CS b u b C la S, 



for 



(0|T{ f^f^JIO), 

7p cs bb ' T c 7a sf 



ST Tr 



+2SflpCS b rC la S: 



2 abc^a'b'c' 
3 fc fc 

+ ST' Ti 



{ 

^cs bb ' T c la s c u c 



+ST Tr 



+2Sf lf3 CSf 1 C<y a S? +2Sf lp CSi b ' 1 C la Sf}, 



(21) 



2 abc a'b'c' i 

3 fc 



1(3 CS bb ' T C la Sf] +S™'Tr \ lf3 CS bb ' T C la Sf 



+2Sf lf5 CS bb C la Sf +2Sfj p CSf C la Sf}, 



(22) 



and for 

(0|T{ Tj?°(x)fjg' o (Q)}\0) i 



ST Tr 



+STi?cs bb c la sf 
+sf 1 ?cs bb ' T c la sf 



2 abc^a'b'c' 
3 fc fc 



{ 



73° Tr 



IpCS, 

T 



+sf lp cs b rc la sz 

+sr'i P cs bb,T c la s<: 



CjaS™ 



+ST Tr 



T 



+STi P cs b rc la S2 
+sf^cs bb ' T c 7a sf}. 



(23) 
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In the above equations, 



Sf(x,0 ] F) = (0\T{q a (x)q b (0)}\0) F , q = u,d,s, (24) 

is the fully interacting quark propagator in the presence of the electromagnetic field. To 
first order in F^ u and m q (assume m u = m d = 0,m s ^ 0), and order x 4 , it is given by 
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We use the convention e 0123 = +1 in this work. The vacuum susceptibilities are defined by 



-|(gg) f/«F Q/3 - -te^F^ \ ^-\ + higher order terms. (25) 



(g<Vg)F = e q x(qq)F pu , 
(q9cG^q) F = e g K(qq)F pu , (26) 
{qgctnv P \G pX l$q)F = ie q i{qq)F pv . 

Note that x has the dimension of GeV -2 , while k and £ are dimensionless. 

The calculation proceeds by substituting the quark propagator into the master formulae, 
keeping terms to first order in the external field and in the strange quark mass. Terms up 
to dimension 8 are considered. The various combinations can be represented by diagrams. 
Fig. [I] shows the basic diagrams considered for the decuplet baryon magnetic moments. 
Fig. |2] shows the diagrams considered for the strange quark mass corrections. Note that 
each diagram is only generic. All possible color permutations are understood. Numerous 
tensor structures emerge from the calculations. Upon ordering the gamma matrices in the 
same order as in the phenomenological side, 18 invariant functions are obtained at the 
corresponding tensor structures. By equating them with those in Eq. (|T6|), QCD sum rules 
are constructed. These invariant functions can be classified by the chirality of the vacuum 
condensates they contain. Eight of them , denoted by WEj, involve only dimension-even 
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condensates, thus we call the corresponding sum rules chiral-even. The other eight, denoted 
by WOj, involve only dimension-odd condensates, and we call the corresponding sum rules 
chiral-odd. Note that previous works such as Refs. |2],f4| use the chirality of the tensor 
structures to refer to the sum rules. The two are opposite. 

To keep the presentation smooth, the complete set of sum rules (a total of 160 for the 
decuplet family) obtained in this work are given in the Appendix in a highly condensed 
form. As it turns out, the validity of a particular sum rule depends on the input parameter 
set. Sum rules that are valid for one set may become invalid for another, and vice versa. For 
this reason, it is useful to present all of the sum rules. Another benefit is that it provides 
a basis for other authors to check the calculation. Sufficient detail is given in this work for 
that purpose. 

The various symbols in the sum rules are explained in the following. The condensate 
parameters are denoted by 

a = -(2tt) 2 (uu), b = (g 2 c G 2 ), (ug c a ■ Gu) = -m 2 (uu). (27) 
The re-scaled current coupling 

\ B = (2n) 2 \ B . (28) 

The quark charge factors e q are given in units of electric charge 

e u = 2/3, e d = -l/3, = -1/3. (29) 

Note that we choose to keep the quark charge factors explicit in the sum rules. The advantage 
is that it can facilitate the study of quark effective magnetic moments. The parameters / 
and (f> account for the flavor symmetry breaking of the strange quark in the condensates and 
susceptibilities: 

= (ss)_ = (sg c a ■ Gs) = = k, = & _ 
{uu) {ug c a ■ Gu) x K £ 

The four-quark condensate is parameterized by the factorization approximation 

(uuuu) = k v {uu) 2 , (31) 

and we will investigate its possible violation via the parameter k v . The anomalous dimension 
corrections of the currents and various operators are taken into account in the leading 
logarithmic approximation via the factor 

(32) 

where \x = 500 MeV is the renormalization scale and Aqcd is the QCD scale parameter. As 
usual, the excited state contributions are modeled using terms on the OPE side surviving 
M 2 — > oo under the assumption of duality, and are represented by the factors 



a s (n 2 ) 
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rin(M 2 /A| CD )- 


7 


k(M 2 )_ 









9 



X 



E n (x) = l-e- x Yl 



x = w 2 B /M 2 B 



(33) 



n 



where wb is an effective continuum threshold. Note that wb is in principle different for 
different sum rules and we will treat it as a free parameter in the the analysis. 

The coefficients for other members of the decuplet family can be obtained by appropriate 
replacements of quark contents. They are: 

1. for A ++ , replace s quark by u quark in Q~ , 

2. A + : replace s quark by d quark in 

3. for A : replace s quark by d quark in 

4. for A~, replace s quark by d quark in Q~ , 

5. for replace u quark by d quark in 

6. for replace u quark by d quark in S*°. 

Here the conversions between u and d quarks are achieved by simply switching their charge 
factors e u and e^. The conversions from s quark to u or d quarks involve setting m s = 0, 
/ = <fi — 1, in addition to the switching of charge factors. 

Furthermore, in the course of collecting the coefficients for the four selected members 
, we discovered some relations among them that allow one to write down 
one set of q staring from another. The relations are given as follows. 

1. From £* + to simply replace every occurrence of e u by (e u + 

2. From H*° to Q~ involves converting the u quark to s quark. This is achieved by 
collapsing each coefficient into a single term that has the maximum number of e s , f, 
(J) in that coefficient. The numerical factor of it is the sum of the numerical factors 
in front each of the terms in the coefficient. For example, (2e s + e u ) goes to 3e s , 
(2e s f(p + e u ) goes to 3e s f(p, (2e s f - 3e s - 3e u f + e u ) goes to -3e s f, etc.. 

These relations were also used as consistency checks of the calculation. 

From the above discussions, we see that it is possible to write down the coefficients for all 
other members of the decuplet family starting just from those for £* + and In the sum 
rule from WEi in Eq. (|A1|) , the complete sets of q are given for the four selected members 
S*° and Q~. They are intended as examples for the reader to get familiar with 
the relations. The rest of the sum rules are presented with q only given for E* + and S*°. 

Finally, let us point out some exact relations among the OPE sides of the sum rules: 



OPE A+ = -OPE A++ , 



(34) 
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OPE A o = 0, 



(35) 



OPE A - = -OPE A+ . (36) 

OPE s »o = l - (OPE E *+ + OPE E *-). (37) 

These results are consequences of symmetries in the correlation functions. As an example, 
let us examine Eq. ([34]). For a given diagram, the master formula for A ++ can be written as 
2{e u C\ + e u C2) = |(Ci + C2) where G\ has the trace dependence, while C2 not. On the other 
hand, the master formula for A + can be written as |[(2e u +ed)Ci + (2e u +ed)C2] = | (C1-I-C2), 
hence the factor of 2. The key here is: a) each term is proportional to a quark charge factor; 
b) SU(2) flavor symmetry in u and d quarks; c) it is the same C\ and C2 that appear in both 
cases. The argument can be generalized to any diagrams, only with C\ and C2 different from 
diagram to diagram. Thus the factor of 2 will survive, regardless of the number of diagrams 
considered. One can argue for the rest of the relations by the same token. The above results 
have been explicitly verified using the calculated coefficients in the sum rules. They also 
provided a set of highly non-trivial checks of the calculation. A number of hard-to-detect 
errors have been eliminated this way. 

Now let us consider the phenomenological side of Eq. ([34]). Since the continuum is 
modeled using terms on the OPE side, the continuum contributions also differ by a factor of 
2. Assuming the transitions, which are modeled by a constant, also differ by a factor of 2, 
then Eq. (p4]) can be extended to the magnetic moments. This assumption was confirmed by 
numerical analysis. The same is true for Eq. ( |3"5] ) and Eq. fl3"6|). The situation for Eq. ( |3"TD 
is a little different. The convergence properties may change when two OPE series are added 
up. Numerical analysis confirmed that fewer sum rules are valid for £*° than for S* + and 



III. MONTE-CARLO ANALYSIS 



To analyze the sum rules, we use a Monte-Carlo based procedure recently developed 
in Ref. ||29|| . The basic steps are as follows. First, the uncertainties in the QCD input 
parameters are assigned. Then, randomly-selected, Gaussianly-distributed sets for these 
uncertainties are generated, from which an uncertainty distribution in the OPE, o 2 OPE {M^) 
where Mj are evenly distributed points in the desired Borel window, can be constructed. 
Next, a x 2 minimization is applied to the sum rule by adjusting the phenomenological fit 
parameters. Note that the uncertainties in the OPE are not uniform throughout the Borel 
window. They are larger at the lower end where uncertainties in the higher-dimensional 
condensates dominate. Thus, it is crucial that the appropriate weight is used in the calcu- 
lation of x 2 - F° r the OPE obtained from the k'th set of QCD parameters, the x 2 P er degree 
of freedom is 

xl 1 ^ [nrW) - nr e "(M j; A fc ,m fc , w k )f 
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where n p is the number of phenomenological search parameters, and n Phen denotes the phe- 
nomenological representation. In practice, ng=51 points were used along the Borel axis. 
The procedure is repeated for many QCD parameter sets, resulting in distributions for 
phenomenological fit parameters, from which errors are derived. Usually, 200 such config- 
urations are sufficient for getting stable results. We generally select 1000 sets which help 
resolve more subtle correlations among the QCD parameters and the phenomenological fit 
parameters. 

The Borel window over which the two sides of a sum rule are matched is determined 
by the following two criteria. First, OPE convergence: the highest-dimension-operators 
contribute no more than 10% to the QCD side. Second, ground-state dominance: excited 
state contributions should not exceed more than 50% of the phenomenological side. The 
first criterion effectively establishes a lower limit, the second an upper limit. Those sum 
rules which do not have a Borel window under these criteria are considered invalid. 



A. QCD Input Parameters 



The QCD input parameters and their uncertainty assignments are given as follows. The 
condensates are taken as 

a = 0.52 ± 0.05 GeV 3 , b = 1.2 ± 0.6 GeV 4 , m\ = 0.72 ± 0.08 GeV 2 . (39) 

For the factorization violation parameter, we use 

k v = 2±1 and 1 < k v < 4. (40) 

The QCD scale parameter is restricted to Aq CT) =0.15±0.04 GeV. The vacuum susceptibili- 
ties have been estimated in studies of nucleon magnetic moments |2|-|J , but the values vary 
in a wide range depending on the method used. Here we take some median values with 50% 
uncertainties: 

X = -6.0 ± 3.0 GeV- 2 and GeV- 2 < x < -10 GeV' 2 , (41) 

and 

k = 0.75 ± 0.38, £=-1.5 ±0.75. (42) 

Note that x is almost an order of magnitude larger than k and £, and is the most important 
of the three. The strange quark parameters are placed at [|5|JT~3[ 



m. 



0.15 ± 0.02 GeV, f = 0.83 ± 0.05, = 0.60 ±0.05. (43) 



These uncertainties are assigned conservatively and in accord with the state-of-the-art in 
the literature. While some may argue that some values are better known, others may find 
that the errors are underestimated. In any event, one will learn how the uncertainties in the 
QCD parameters are mapped into uncertainties in the phenomenological fit parameters. In 
the numerical analysis below, we will also examine how the spectral parameters depend on 
different uncertainty assignments in these input parameters. 
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B. Search Procedure 



To extract the magnetic moments, a two-stage fit was performed. First, the correspond- 
ing chiral-odd mass sum rule, as obtained previously in Ref. |U|, was fitted to get the mass 



Mb, the coupling X 2 B and the continuum threshold w\. Then, Mb and were used in the 
magnetic moment sum rule for a three-parameter fit: the transition strength A, the contin- 
uum threshold W2, and the magnetic moment fi B . Note that W\ and W2 are not necessarily 
the same. We impose a physical constraint on both w± and W2 requiring that they are larger 
than the mass, and discard QCD parameter sets that do not satisfy this condition. In the 
actual analysis of the sum rules, however, we found that a full search was not always suc- 
cessful. In such cases, the search algorithm consistently returned W2 either zero or smaller 
than Mb- This signals insufficient information in the OPE to completely resolve the spec- 
tral parameters. To proceed, we fixed W2 at wi, which is a commonly-adopted choice in the 
literature, and searched for A and fi B . The two-stage fit incorporates the uncertainties from 
the two-point functions in a correlated fashion into the three-point functions, and represents 
a more realistic scenario. 

To illustrate how well a sum rule works, we first cast it into the subtracted form, 

U s = ~X 2 B^e' M B/ M \ (44) 



then plot the logarithm of the absolute value of the two sides against the inverse of M 2 . In 
this way, the right-hand side will appear as a straight line whose slope is — M B and whose 
intercept with the y-axis gives some measure of the coupling strength and the magnetic 
moment. The linearity (or deviation from it) of the left-hand side gives an indication of 
OPE convergence, and information on the continuum model and the transitions. The two 
sides are expected to match for a good sum rule. This way of matching the sum rules is 
similar to looking for a 'plateau' as a function of Borel mass in the conventional analysis, but 
has the advantage of not restricting the analysis regime in Borel space to the valid regimes 
common to both two-point and three-point correlation functions. 



IV. RESULTS AND DISCUSSIONS 



We have analyzed all of the sum rules for the entire decuplet family. We confirmed the 
three relations among magnetic moments as extended from Eqs. ( |3"4"D to (^). So we will 
only present results for seven members. Valid sum rules were identified using the criteria 
discussed earlier. The results are given in three tables: Tables | to [TIT]. The corresponding 
overlap plots are given in seven figures: Figs. ^| to |9|. These plots show how well a sum 
rule performs in the entire Borel region. Such information is absent in the tables. From the 
results, the following observations are in order. 

In general, more chiral-even sum rules are valid than chiral-odd ones. This is consistent 
with previous findings for the octet baryon magnetic moments. It was argued in Ref. || that 
the interval of dimensions (not counting the dimension of F^) in the chiral-even sum rules 
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(0 to 8) is larger than that in the chiral-odd sum rules (1 to 7). Indeed, more chiral-even 
sum rules (WE 2 , WE 4 , WE 5 , WE 6 , WE 8 , WE 9 ) have power corrections up to 1/M 4 , than 
chiral-odd ones (W0 2 , W0 4 , W0 8 , W0 9 ). Because of the additional terms in the OPE 
series, these sum rules are expected to be more reliable than the other sum rules. The 
situation here is almost opposite to that for the two-point functions flI3|. It was pointed out 



in Ref. |30[ that chiral-odd sum rules are more reliable than chiral-even sum rules for baryon 
two-point functions. The reason could be traced to the fact that even and odd parity excited 
states contribute with different signs. In the three-point functions, however, the statement is 
no longer valid due to the appearance of transitions and vacuum susceptibilities. Therefore, 
caution should be used when applying the chirality argument to determine the reliability of 
a sum rule in three-point functions. In addition, numerical analysis showed that the sum 
rules from WEi, WE 3 , WOi are valid for the standard input parameter set, despite the 
absence of 1/M 4 terms. We have varied the central values of the input parameters and 
discovered that sum rules that were valid for one set of input parameters became invalid for 
another, and vice versa. Thus the situation with three-point functions is more complicated. 
Our experience is that each sum rule should be examined individually in order to find out 
its reliability. 

It turns out that for most of the valid sum rules, a full search was unsuccessful, except 
for three sum rules: WE5, WEg and WEg for Of the three, only WEg returned a 

continuum threshold with reasonable error. The other two returned it with large errors. 
The large errors indicate that the sum rules are not very stable: they contain barely enough 
information to completely resolve the spectral parameters.. The important point is that the 
results with the continuum threshold searched or not are almost the same. This suggests 
that fixing it to that of the corresponding two-point function seems a good approximation. 

It is gratifying to observe that the valid sum rules for most members give consistent 
predictions for the magnetic moments in terms of the sign, except for E*° and S*° whose 
magnitudes are small. The magnitudes for the magnetic moments are consistent within 
errors for the most part, with only a few exceptions. The performances of the sum rules are 
quite different within each member. This is best displayed in the overlap plots. In some 
sum rules, the overlap is poor, as evidenced by the deviation from linearity (dotted lines). 
It signals poor OPE convergence in these sum rules. As expected, the deviation is more 
severe in the lower end of the Borel region where nonperturbative physics dominates. These 
sum rules will more likely suffer from uncertainties associated with the selection of the Borel 
window. As a result, the spectral parameters extracted from them are less reliable. One way 
to alleviate the problem is to increase the lower end of the Borel window to values where 
the overlap is good, even to extend the upper end to ensure the existence of a window. This 
was not attempted in this work because we feel the results obtained this way are somewhat 
misleading. The reason is that the sum rules in these windows will be dominated mostly 
by perturbative physics. It is common knowledge that if one goes deep enough into the 
Borel space, one can always find a match in a QCD sum rule. But such practice is against 
the philosophy of the QCD sum rule approach, which relies upon the power corrections to 
resolve the spectral properties. Therefore, some standard is necessary to emphasize such 
physics, and we feel the 10%-50% criteria adopted here are a reasonable choice. 
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Based on the quality of the overlap, the broadness of the Borel window and its reach 
into the lower end, the size of the continuum contribution, and the standard QCD input 
parameter set, we designate one sum rule for each member as the most favorable. They are 
WE 5 for A++, WE 5 for £*+, W0 8 for WE 5 for E*~ W0 8 for WE 5 for E*~, and 
WE5 for Q~ . The selection is undoubtedly subjective. The reader may find a different set 
that have equal or comparable performance. We want to stress that such a selection depends 
on the QCD input parameters. It is possible that the ones selected here become invalid for 
a different set of input parameters, in which case a new set should be selected. 

Relatively large errors were found in the valid sum rules using the standard QCD input 
parameter set: from 50% to 100% in the magnetic moments. But in most cases, the sign 
and order of magnitudes are unambiguously predicted when compared to the measured 
values. The situation is similar to a previous finding on g& fl3l| . To gain some idea on 
how the uncertainties depend on the input, we also analyzed the sum rules by adjusting the 
error estimates individually. We found large sensitivities to the quark condensate magnetic 
susceptibility x- I* 1 fact, most of the errors came from the uncertainties in \. We also 
tried with reduced error estimates on all the QCD input parameters: 10% relative errors 
uniformly. It leads to about 30% accuracy on the magnetic moments in the favorable sum 
rules. Further improvement of the accuracy by reducing the errors in the input is beyond the 
capability of these sum rules as the \ 2 '/Ndf becomes unacceptably large, signaling internal 
inconsistency of the sum rules. For that purpose, one would have to resort to finding sum 
rules that have better convergence properties and depend less critically on the poorly-known 
X- 

To get a different perspective on how the spectral parameters depend on the input pa- 
rameters, we study correlations among the parameters by way of scatter plots. In the 
Monte-Carlo analysis, all the parameters are correlated. Therefore, one can study the corre- 
lations between any two parameters by looking at their scatter plots. Such plots are useful in 
revealing how a particular sum rule resolves the spectral properties. We have examined nu- 
merous such plots. Here we focus on the favorable sum rules as selected earlier. To conserve 
space, we only give two examples. Fig. ITU] shows the scatter plot for correlations between 



Q~ magnetic moment and the QCD input parameters for the sum rule from WE5. Fig. |Tl 
shows a similar plot for S*° and the sum rule from WOs- Perhaps the most interesting 
feature is the strong correlations with x m both sum rules. This is the reason for the large 
sensitivities to this parameter as alluded to earlier. Precise determination of \ is crucial for 
keeping the uncertainties in the spectral parameters under control. Other charged members 
(all use sum rules from WE 5 ) display qualitatively the same patterns for parameters other 
than x an d the factorization violation parameter k v . For x, positively-charged members 
(A ++ and show negative correlations. The opposite is true for negatively-charged 

members S*~ and S*~): they show positive correlations with x- The patterns for k v 
essentially follow those for x, although the correlations are weaker. The correlation patterns 
for S*° are qualitatively the same as those for 

Table [TV] shows a comparison of the magnetic moments from various calculations and 
existing experimental data. The results with 10% errors from the QCD sum rule method 
are used in the comparison. Note that the central values are slightly different from those in 
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Tables | to [Hj where conservative uncertainties were used. The reason is that the resultant 
distributions vary with input errors and are not Gaussian in this case. In such event the 
median and the average of the asymmetric errors are quoted. The QCDSR results are 
consistent with data, although the central value for Q~ is slightly underestimated. We 
would like to point out that it is possible to reproduce the central value for Q~ (using it as 
input) by fine-tuning of the susceptibility \ alone, given the sensitivity to this parameter and 
the large freedom at the present time on its value. However, we feel that such an attempt is 
not very meaningful given the accuracy of the method. A more meaningful practice would 
be to re-analyze the octet baryon magnetic moments by the same method as employed here, 
and obtain a best fit on the the susceptibilities using their accurately measured values, then 
use them to predict the decuplet magnetic moments. It would yield valuable information on 
these important quantities and on the consistency of the approach. From the table, it is fair 
to say that the QCDSR approach is at least competitive with other calculations. The results 
came about from a rather different perspective: the nonperturbative structure of the QCD 
vacuum. The results from various calculations roughly agree, except for the charge- neutral 
resonances A , and ~*° for which both the sign and the magnitude vary. It would be 
helpful to have experimental information on the other members of the decuplet, although 
such measurements appear difficult. 



V. CONCLUSION 



It has been demonstrated in this work that the magnetic moments of decuplet baryons 
can be successfully computed in the QCD sum rule approach. A complete set of QCD sum 
rules are derived using the external field technique. They are analyzed extensively with 
a comprehensive Monte-Carlo based procedure which , in our opinion, provides the most 
realistic estimates of the uncertainties present in the approach. 

Valid sum rules are identified using criteria established by OPE convergence and ground- 
state dominance. For each member, usually several sum rules are valid, but not all of 
them perform equally well. This was best displayed by the overlap plots. Some have large 
deviations in the lower end of the Borel window, signaling insufficient convergence in the 
OPE. These sum rules are less reliable. Based on overall performance, a favorable sum 
rule was selected for each member. They are WE5 for charged members, WOg for charge- 
neutral members. We also found the following relations between the magnetic moments: 
A*a+ — |A i A++ ) Ha — 0, and /i A - = — fj, A +, and approximately /i E .o = |(yU E .<- + + /i E .-). . 

Using conservative estimates of the QCD input parameters, the uncertainties in the ex- 
tracted magnetic moments are found relatively large as compared to the two-point functions. 
We found that the results are sensitive to the quark condensate magnetic susceptibility x- In 
fact, most of the uncertainties could be attributed to x- Better estimate of this parameter is 
clearly needed. By varying the uncertainty estimates in the input parameters, we found that 
a 30% accuracy can be achieved with the designated sum rules if the QCD input parameters 
could be determined to the 10% accuracy level. 
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APPENDIX A: QCD SUM RULES FOR MAGNETIC MOMENTS OF 

DECUPLET BARYONS 

Here we give the complete set of QCD sum rules derived in this work. Foe each member, 
there are 18 sum rules, 9 chiral-even, 9 chiral-odd. It turns out that the sum rules from WE 6 
and WE 7 are degenerate, so are those from W0 4 and W0 5 . So the number of independent 
sum rules is 16 for each member. The total number for the entire decuplet family is 160. 
They are given in the following in a highly compact form. The explanation on how to obtain 
a sum rule for a particular member is discussed in the main text. 

The sum rule from WE^ 

c 1 L 4 / 27 E 1 M 4 + c 2 m sX aL- 12 ' 27 E M 2 + c 3 bL^ 27 + c iX a 2 L 12 / 27 + (c 5 + c 6 )m s aL 4 / 27 
+ (c 7 + c 8 )a 2 L 2 ^ 27 ^ + c 9X m 2 a 2 L- 2 / 27 ^ + c w m s m 2 aL^/ 27 ± 

= \^(w* + A ) e ~ MllM2 ' (A1) 



where the coefficients for S* + are: 

ci = |(e s + 2e u ), c 2 = ^(e s /0 + 2e u ), 

c 3 = ^ (e s + 2e u ) , c 4 = =± (e s f<p + e u f + e u ) , 

c 5 = ±(-2eJ + 9e s + 9e u f + 5e u ), c 6 = =±{e 8 f<f> + 2e M )(7/t + 0, (A2) 

c 7 = M~ e *f + 3e * + 5e «/ _ e u )««» c 8 = (e a /0 + e u f + e u )(7n + £), 

c 9 = 2Tfi(e s /0 + e u f + e„), c w = (e s + e u f + e u ), 



for S 



ci = |(e s + e u + e d ), c 2 = -^(e s f<f> + e u + e d ), 

c 3 = ^(e s + e u + e d ), c 4 = ^{e s f<p+ (e u + e d )(f + l)/2), 

c 5 = ±(-2e s f + 9e s + (e u + e d )(9f + 5)/2), c 6 = =±{e a f<f> + e u + e d )(7 K + £), 

c 7 = U-eJ + 3e s + (e u + e d )(5f - 1)/2)k v , c 8 = ^(e 8 f<f> + (e u + e d )(f + 1)/2)(7k + 0, 

c 9 = 2l6( e ^ + ( e « + e <i)(f + i)/ 2 )' c io = rf (e s + (e„ + e d )(/ + l)/2), 

(A3) 

for S*°: 
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for Q~ 



ci = |(2e s + e u ), c 2 = ^f(2e s /0 + e u ), 

c 3 = ^-(2e s + e u ), c 4 = ^r/(e s /0 + e s + e n ), 

c 5 = ^(5e s / + 9e s + 9e u /-2e u ), c 6 = ^(2e s /0 + e u )(7« + 0, (A4) 

c 7 = j? f(-e s f + 5e s + 3e M / - e M )/t„, c 8 = ^f(e s f(j) + e s (f) + e u )(7/c + 0, 

c 9 = ?kf(. e sf<P + e s + e u ), cio = (e s / + e u / + e s ), 



ci = |e s , c 2 = =j-e s f(f), 

c 3 = ii e s, c 4 = -y-e s / 2 0, 

C5 = |e s /, c 6 = =±e s f(f){7K + 0, (A5) 

C7 = ^e s f 2 K v , c 8 = ^e s f 2 0(7/t + £), 

c 9 = ^e s f 2 (j), cio = if e a /. 



The sum rule from WE 2 : 



Cl L^ 27 E M 2 + c 2 m sX aL-^ 27 + c 3 bL^ 27 ^ + c 4 m s aL^ + c 5 m s m 2 aL^ 27 ^ 



9 B I Ml M 2 



% A + ^ e"^, (A6) 



where the coefficients for S* + are: 

ci = ^(e s + 2e u ), c 2 = ^r(e s /0 + 2e M ), c 3 = -^(e s + 2e„), c 4 = |(e s + e u / + e u ), 

for 

ci = ^-(2e s + e tt ), c 2 = =p(2e s /0 + e u ), c 3 = ^-(2e s + e u ), c 4 = |(ej + e u / + e s ), 
c s = Ts( e sf + e uf + e s ). 

The sum rule from WE 3 : 

Cl L A l 21 E x M A + c 2 m sX aL- l2 l 27 E M 2 + c 3 6L 4 / 27 + c 4X a 2 L 12 / 27 + (c 5 + c 6 )m s aL A / 27 
+ (c 7 + c 8 )a 2 L 28 / 27 ^ + c, X mla 2 L- 2 l 27 ±- 2 + c 10 m s m 2 aL- 10 / 27 _L 

= lJ^(^ + A ) e " M|/M2 ' (A9) 

where the coefficients for E* + are: 

ci = ^(e s + 2e u ), c 2 = ^(e s f(f) + 2e u ), 

c 3 = ^g(e s + 2e M ), c 4 = jg(e s f<p + e u / s + e u ), 

c 5 = ^(2eJ-3e s -3e n / + e u ), c 6 = i(e + 2e u )(4« + 0, (A10) 

c 7 = ^( e s/ - 3e s /2 - 2e u / + e u )/c„ c 8 = ^{e s f4) + e «/ + e„)(4« + 0, 

eg = ^(e s f(f) + e u f + e u ), c l0 = ^(e s + e u f + e u ), 
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for S 



c i = 24 ( 2e « + e «)' c 2 = 2i(2e s /0 + e n ), 

c 3 = g|g(2e s + e u ), c 4 = ^f{e s f(f) + e s (p + e u ), 

c 5 = ^(e s / - 3e s - 3e u f + 2e u ), c 6 = ^(2e s /0 + e M )(4/t + (All) 

c 7 = hf( e sf - 2e s - 3e M //2 + e u )/c„ c 8 = ±f{e 8 f<j> + e s + e u )(4/c + f), 

The sum rule from WE4: 

CiL 4/27 £!M 4 + c 2 m s xa^~ 12/27 ^ M 2 + c 3 6L 4/27 + (c 4 + c 5 )m s aL 4/27 
+ (c 6 + c0a 2 L 28 / 27 ^ + c 8 m 2 a 2 L 44 / 27 I L 



18 B VM 2 J 



c 5 

c 7 = TM( e */0 + e uf + e«)(4/t + 0, c 8 = ^(e s f + e u f + e u ), 



for S 



(A12) 



where the coefficients for £ are: 

ci = ^(e s + 2e u ), c 2 = ^(e s f(p + 2e u ), 

C3 = i^(e s + 2e M ), c 4 = j^(e s f + Qe s + Qe u f + 8e u ), . 

^(e s /0 + 2e tt )(10« + 0, Q> = £(e a / + 3e s /2 + 4e u / + e„K, 1 JJ 



ci = ^(2e s + e u ), c 2 = ^|(2e s /0 + e u ), 

c 3 = g^(2e s + e u ), c 4 = yg(8e s / + Qe s + Qe u f + e u ), , . 

c 5 = ^(2e 8 /0 + e u )(lOK + O, c 6 = £/(e J + 4e s + 3e u //2 + e u K, 1 J 

c 7 = mfi^f^ + e ^ + e «)( 4/t + 0» c s = ^f( e sf + e s + e u ). 

The sum rule from WE5: 

c^E.M 4 + c 2 m sX aL- 12 / 27 £ M 2 + c 3 bL^ 27 + (c 4 + c 5 )m s aL^ 27 + c 6X a 2 L 12 ' 27 

+ (c 7 + c 8 )a 2 L 28 / 27 ^ + C9Xm 2 a 2 L- 2 / 27 ^ + c^m 2 ^- 10 / 27 ^ + c u m 2 a 2 L^ 27 ^ 



18 B VM 2 / 
where the coefficients for S* + are: 



(A15) 



c i = 2l( e s + 2e u ), c 2 = 1 |(e s /0 + 2e u ), 

c 3 = iii( e s + 2e M ), c 4 = |(e 8 / + e s + e„/ + 3e tt ), 

C5 = ^{e s f<P + 2e «)( 2/t + 0> c 6 = 9( e J0 + e «/ + e «)> / Al6 x 

c 7 = ^(e a / + e u / + e u )n v , c 8 = j^(e s f(f) + e u f + e u )(4n + f), 

eg = ^(e s f(f) + e u f + e u ), c 10 = ^(e s + e u f + e u ), 
°n = ^( e sf + e u f + e„), 
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for S 



ci = 2j(2e s + e M ), c 2 = 

c 3 = Tii( 2e s + e «)> C 4 = 

c 5 = ±(2e s f(fr + e u )(2n + f ), c 6 = 

c 7 = §ff(e s f + e s + e u )n v , c 8 = 

c 9 = 2^f(e s f(f) + e s (f) + e u ), c 10 
c n = e^/(e s / + e s + e u ). 



«7 5 



|(3e s / + e s + e u f + e u ), 



6 

T^8/(e s /0 + e s< 



+ e tt )(4« + 0, 



(A17) 



j_ 

24 



e s / + e u f + e 8 ), 



The sum rule from WE 6 : 



Cl m sX aL- l2 l 27 E M 2 + (c 2 + c 3 )m s aL 4 / 27 + (c 4 + Cg)^ 28 / 27 --^ 
+c 6 m s m 2 aL- 10 / 27 ^ + c 7 m 2 a 2 L^ 27 ^ 

= ^ I (^ + A ) 6 " M|/M2 ' (A18) 



where the coefficients for S* + are: 

ci = j L (e s f(j) + 2e u ), c 2 = j^(e s f + 2e u ), 

c 3 = ^(e s /0 + 2e u )(2/c - £), c 4 = ^(e a / + e„/ + e u )/t t 

c 5 = i5s( e «/0 + e «/ + e «)( 2/t - 0; c 6 = ^( e « + e «/ + e «)' 



for S*°: 



ci = ^(2e s /0 + e M ), c 2 = j^(2e s f + e u ), 

c 3 = ^(2e s /0 + e M )(2/t - £), c 4 = ^f(e s f + e s + e u )n v , 

°7 = ^f(e s f + e s + e u ). 



(A19) 



(A20) 



The sum rule from WE 7 is identical to that from WE 6 after multiplying an overall sign 
on both sides. 

The sum rule from WE 8 : 

crnisXaL- 12121 E M 2 + c 2 bL 4/27 + c 3X a 2 L 12/27 + (c 4 + c 5 )m s aL 4/2J 

+ (c 6 + c 7 )a 2 L 28 / 27 ^ + C8Xm 2 a2jL - 2/27 _l_ + C9msm 2 aL -io/ 27 _l_ + Clom 2 a2jL i4/ 27 _l_ 

_-2~ 

~ 3 VM^ 

where the coefficients for S* + are: 

ci = i(e s /0 + 2e u ), c 2 = itAe s + 2e u 



(^ + Ae- M ° /M \ (A21) 



3 \ sj r i — «/! 144 \ 

c 3 = y (eJ0 + e tt / + e M ), c 4 = =±(5e a f + 12e s + 12e u f + 22e u ), 

c 5 = |(e a /0 + 2e tt )(14K-O, c 6 = ^(8e s / + 6e s + 20e u / + 8e u )n v , (A22) 

c 7 = ^(e s f(f) + e u f + e u )(4/c + f ), c 8 = YSg( e J0 + e «/ + e «)> 

eg = ^(e s + e u f + e u ), c w = m (e s f + e u f + e u ), 
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for 



ci = l(2e s f<f> + e u ), c 2 = ■ 

c 3 = ^/(e s /0 + e s + e u ), c 4 = ■ 

c 5 = ^(2e s /0 + e u )(14/c-O, Q> = 

c 7 = uf(e s f(f) + e s (f) + e u )(4/c + f), c 8 = • 

c 9 = ^(e s / + e u / + e s ), c 10 = 



144 (^Cs + C u ), 

^(22eJ + 12e s +12e u / + 5e n ), 
^/(8eJ + 20e s + 6e u / + 8e„K, (A23) 
lM/(e s /0 + e s + e u ), 
: iHg/(eJ + e s + e u ). 



The sum rule from WE 9 has the same form as that from WE 8 , only with different cf 

ci = (e s f(j) + 2e u ), c 2 = j^(e s + 2e u ), 

c 3 = |(eJ0 + e u / + e u ), c 4 = |(4eJ + 12e s + 12e u f + 20e M ), 

c 5 = ^(eJ0 + 2e u )(-8n + £), c 6 = |(2eJ + 2e s + 6e u / + 2e u )n v , (A24) 

C7 = ^r(e s /0 + e u / + e u )/c, c 8 = ^(e 8 f<f> + e u / + e„), 

c 9 = ^(e s + e„/ + e u ), ci = ^(e s / + e u / + e u ), 



ci = (2e s /0 + e u ), c 2 = ^(2e s + e u ), 

c 3 = |/(eJ0 + e s + e u ), c 4 = |(20eJ + 12e s + 12e u f + 4e„), 

c 5 = ^(2eJ0 + e u )(-8K + £), c 6 = \f{2e s f + 6e s + 2e u f + 2e„K, (A25) 

C7 = ir/(e s /0 + e s + eu)«, c 8 = j£f(e a f<f> + e s <p + e„), 

c 9 = =±(e s f + e u f + e s ), c 10 = f^f(e s f + e s + e u ). 



The sum rule from WOi: 

cixaEiM 4 + c 2 m s L m/27 E X M 4 + (c 3 + c 4 )aL 16/27 E M 2 + c 5 m 2 aL 2/27 

+CQX a b + c-jm s xa> A 
'H B M B 



CsabL^jL + (c 9 + c 10 )m s a 2 L 16 / 27 ^ 



2 B 



M 2 



M|/M 2 



(A26) 



where the coefficients for X are: 



Cl = 


^(eJ0 + 2e u ), 


c 2 = • 


c 3 = 


~2f( e sf ~\~ 2e u ), 


c 4 = • 


c 5 = 


12 Cm,/ C«) ? 


c 6 = ■ 


c 7 = 


y(e s /0 + e u / + e tt ), 


c 8 = • 


c 9 = 


2j(-2e.J + 3e s + 4e u / - 2e u )/c t; , 


ClO = 



j2(eJ0 + 2e u )(-8 
^(e s /0 + 2e u ), 



432 1 



(A27) 









Cl = 


^(2eJ0 + e u ), 


c 2 = 


c 3 = 


=f(2e s / + e u ), 


c 4 = ■ 


c 5 = 


^(e s / + e u / + e s ), 


c 6 = 


c 7 = 


y/(e s /0 + e s + e u ), 


c 8 = ■ 


c 9 = 


^/(-2e s / + 4e s + 3e u /- 


2e u )K v , cio = 



g (2c s "I - c u ) , 

^(2e i /0 + e tt )(-8K + 7O, 
^(2e s /0 + e u ), 

2ig( e s/ + Cm/ + C s ), 

^f(e s f(j) + e s + e u )(8/c + 



(A28) 
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The sum rule from WO2: 



c^aEoM 2 + (c 2 + c 3 )aL w/27 + c 4 m sX a 2 + c^xah^— + (c 6 + c 7 )m s a 2 L w/27 



[J. 1 



Mb M 2 



A) e~ M H M \ 



where the coefficients for S* + are: 



ci = |(e s /0 + 2e M ), c 2 
c 3 = ^(eJ0 + 2e u )(5/t + 2O, c 4 

c s = 2i6( es /^ + 2e «)' 
C7 = ir(e s /0 + e u / + e u )/c, 



27 (^s/ ~l~ 2e u ) , 
: f(e s f(j) + e u f + e u ), 



for S 



ci = \{2e s f(j) + e u ), c 2 
c 3 = ^(2eJ0 + e u )(5/t + 2O, c 4 
C 5 = 2l6( 2e s/0 + e u ), c 6 



^(2e s / + e„) 



—f(e s f(j) + e s (f) + e u ), 
~irf(, e sf + e s + e u )/t„, 



The sum rule from W0 3 : 

cixa^iM 4 + (c 2 + c 3 )aL ie/27 E M 2 + c 4 m 2 aL 2/27 + c^xab + c 6 m s xa 2 
+c 7 afeL 16 / 27 ^ + (c 8 + c 9 ) ms a 2 L 16 / 27 -i 



7~ 2 //i s M 



M 2 



B + A) e- M l' M2 



where the coefficients for S* + are: 



Cl = 


^(eJ0 + 2e M ), 


c 2 = 


c 3 = 


^(ej^ + 2e tt )(2«-ll0, 


c 4 = 


c 5 = 


i72|(e s /0 + 2e u ), 


c 6 = 


c 7 = 


1728 

432(e s /0 + e u / + e tt )(12/c + 70, 


c 8 = 


Cg = 





^(7eJ + 6e s + 6e u / + 20e n ), 

1 2 \&s "I - Cm/" "I - C u ), 

^(eJ0 + e M / + e n ), 
"3g"(c s / + 3e s + 7e M / + 6 u )k v , 



for S 



K2e s 



+ e u ) 



c i = h 

c 3 = 2i8(2e s /0 + e u )(2K-llO, 

C 5 = T7l|( 2e s/0+ e n), 

C 7 = T^8( e J + e «/ + e S )> 

c 9 = ^f{e s f<t> + e ^ + e u)(l2K + 70- 



c 2 = - 7 \(20e s f + 6e s + 6e u f + 7e M ), 
c 4 = =^{e s f + e tt / + e s ), 
ce = T2f(e s f<fr + e s (p + e tt ), 



The sum rule from W0 4 : 
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Cl m s L- s/27 E M 2 + (c 2 + c 3 )aL 16 / 27 + c A m sX a 2 + c 5 m 2 aL 2 / 27 — + c %X ab^ 
+o T abLW± i + ( C8 + c 9 )m s a 2 L^ 27 ^ 

= ^b(tJ^ + A^ /M \ (A35) 



18 " \M B M 2 

where the coefficients for S* + are: 

c i = \( e s + 2e u ), c 2 = |(e s + e u / + e tt ), 

c 3 = Tii(e s /0 + 2e M )(12/t + £), c 4 = g(e s /0 + e„/ + e u ), 

c 5 = ij(e s + e tt / + e u ), c 6 = §g(e a /0 + 2e u ), (A36) 

c 7 = 4jb( e * + e «/ + e «)> c s = lk( e s + 2 e u f)K v , 
eg = 432(e s /0 + e u f + e u )(12/c + f), 

for 

ci = |(2e s + e u ), c 2 = |(e a / + e tt / + e s ), 

c 3 = T4i(2e s /0 + e u )(12/t + O, c 4 = gf(e s f(f> + e s (f) + e u ), 

c 5 = T^(e s / + e u f + e s ), c 6 = ^{2e s f(p + e u ), (A37) 

c 7 = ib( e s/ + e «/ + e «)> c s = Ti"/( 2e s + e «/K> 

C 9 = ^fM^ + e ^ + e «)( 12/t + 0- 

The sum rule from W0 5 is identical to that from W0 4 after multiplying an overall sign 
on both sides. 

The sum rule from W0 6 : 

Cixa^iM 4 + c 2 m s L- 8/27 ' E X M^ + (c 3 + c 4 )aL w/27 E M 2 
+c 5 mlaL 2 / 27 + c 6X ab + c 7 abL^ 27 ^ + (c 8 + c 9 )m s a 2 L^ 27 ^ 

= ^(^ + A ) 6 " M|/M2 ' (A38) 

where the coefficients for S* + are: 

ci = ^-(e s /0 + 2e u ), c 2 = =^(e s + 2e u ), 

c 3 = #(lle s / + 18e s + 18e tt / + 40e M ), c 4 = ^(e s /0 + 2e u )(14/c + 23£), 

c 5 = |(e s + e u / + e u ), c 6 = ^(e s f(p + 2e u ), (A39) 

c 7 = ijj|( e « + e «/ + e «)' c § = ^(5e s f + 3e s + lie,,/ + 5e u )K„, 

c 9 = ^( e s/0 + e n/ + e n )(2K - 30, 

for S*°: 

ci = ^-(2e s /0 + e u ), c 2 = ^-(2e s + e u ), 

c 3 = f±(40e s / + 18e s + 18e u / + lle„), c 4 = ^(2e s /0 + e tt )(14/c + 23£), 

c 5 = |(e s / + e u / + e s ), c 6 = ^ ( 2e »/0 + e «)> ( A4 0) 

c 7 = i|( e s/ + e «/ + e «)' c § = 4?f( 5e *f + lle s + 3e «/ + 5e «)^> 

c 9 = iifM<f> + e s (p + e u )(2n - 3f )• 
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The sum rule from WO7 has the same form as that from WC>6 after multiplying an overall 
sign on both sides. They only differ in c 4 and c 9 for 

c 4 = 2ifi(eJ0 + 2e «)( 14/t + 29 ^' c 9 = I^ e ^ + e ^ + e «)(- 4/t + 5 0> (A41) 

for 

c 4 = ^( 2e ^ + e tt )(14« + 290, c 9 = ^/(eJ0 + e s + e„)(-4« + 5£). (A42) 

The sum rule from W0 8 : 

c lX aE Q M 2 + c 2 m s L- 8/27 E M 2 + (c 3 + c 4 )aL 16/27 + c 5 m 2 aL 2/27 ^ 
+c 6 xabj^ + c 7 abL w / 27 ^ + (c 8 + c 9 )m s a 2 L 16 / 27 ^ 

= ( _J^_ + ^ e -^i/^ 2 , (A43) 

3 B \M B M 2 J K ' 

where the coefficients for S* + are: 

ci = ^(e s /0 + 2e u ), c 2 = ^(e^ + 2e u ), 

c 3 = f^(5e s / + 9e s + 9e u f + 19e u ), c 4 = ^{e s f<j> + 2e u )(8/t - £), 

c 5 = |(e s + e u / + e u ), c 6 = ^ (e s /0 + 2e M ), (A44) 

c 7 = 2i|( e s + e «/ + e «)> c § = |( 2e s/ + e s + 4e u / + 2e u )«„, 

c 9 = mi^ffi + e «/ + e «)( 6/t - 0> 

for 

ci = =±(2e a f<f> + e u ), c 2 = =±(2e s + e u ), 

c 3 = ^(19e s / + 9e s + 9e u f + 5e u ), c 4 = ^g(2e s f(f> + e u )(8/t - £), 

c 5 = |(e s / + e u / + e s ), c 6 = ^{2e s f<f) + e u ), (A45) 

c 7 = 2i|( e «/ + e «/ + e «)> c « = lf( 2e sf + 4e s + e u / + 2e u )K„, 

c 9 = mfMJ* + e ^ + e «)( 6/t _ 0- 

The sum rule from WO9 has the same form as that from WOg. They only differ in c 4 
and c 9 for 

C4 = ^(e s f(f) + 2e u )(2K- £), c 9 = -^(e s f(p + e u f + e u )n, (A46) 



for S*°: 



C4 = -^(2e s /0 + e u )(2K-£), c 9 = -^f(e s f(f) + e s (p + e u )n. (A47) 
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TABLE I. Monte-Carlo analysis of the QCD sum rules for the magnetic moment of A ++ and 
Q~. The six columns correspond to, from left to right: the sum rule that has a valid Borel region, 
the Borel region determined by the 10%-50% criteria, the percentage contribution of the excited 
states and transitions to the phenomenological side at the lower end of the Borel region (it increases 
to 50% at the upper end), the continuum threshold, the transition strength, the magnetic moment 
in nuclear magnetons. The uncertainties in each sum rule were obtained from consideration of 1000 



QCD parameter sets. 


Sum Rule 


Region 


Cont 


w 


A 


Ms 




(GeV) 


(%) 


(GeV) 


(GeV" 2 ) 




A++: WEi 


0.70 to 1.53 


1.7 


1.65 


-0.28 ± 0.52 


7.76 ± 2.67 


WE 3 


1.04 to 1.42 


19 


1.65 


0.20 ± 0.20 


3.06 ± 1.14 


WE 4 


0.675 to 1.56 


5 


1.65 


-0.35± 0.37 


3.34 ± 1.44 


WE 5 


0.765 to 1.47 


8.5 


1.65 


0.53± 0.81 


3.56± 3.49 


Sl~: WEi 


0.592 to 1.70 


2.3 


2.30 


-0.12± 0.11 


-2.66 ± 0.88 


WE 2 


0.872 to 1.53 


20 


2.30 


-0.26± 0.20 


-5.31 ± 3.66 


WE 3 


0.885 to 1.68 


8.5 


2.30 


-0.09± 0.04 


-1.24 ± 0.51 


WE 4 


0.60 to 1.72 


2 


2.30 


-0.03± 0.05 


-1.24 ± 0.24 


WE 5 


0.747 to 1.66 


7.4 


2.30 


-0.14± 0.14 


-1.32 ± 1.08 


WE 6 


0.59 to 2.32 


0.86 


2.30 


-0.01 ± 0.02 


-1.14 ± 0.40 


WE 8 


0.69 to 2.60 


3.3 


2.30 


0.03 ± 0.03 


-0.65 ± 1.22 


WOi 


0.663 to 1.26 


12 


2.30 


-0.32± 0.42 


-0.65 ± 1.22 


wo 2 


1.06 to 1.43 


31 


2.30 


-0.62± 0.18 


-4.94 ± 5.58 


wo 4 


0.836 to 2.22 


7.4 


2.30 


-0.03± 0.01 


-0.70 ± 0.24 
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TABLE II. 


Same as Table |, 


but for £*+ 


, S*° and E*~. 


The presence of a second row in a 


specific sum rule indicates that the continuum threshold was 


successfully searched. 




Sum Rule 


Region 


Cont 


w 


A 


Ms 




(GeV) 


(%) 


(GeV) 


(GeV~ 2 ) 




WEi 


0.853 to 1.445 


11 


1.80 


0.28± 0.18 


2.96 ± 1.41 


I » J — « Q 


0.996 to 1.39 


23 


1.80 


0.23± 0.08 


1.49 ± 0.79 


WE 4 


0.622 to 1.61 


1 


1.80 


-0.05± 0.10 


1.74 ± 0.42 


WE^ 


0.715 to 1.45 


10 


1.80 


0.34± 0.35 


1.82 ± 1.94 




0.715 to 1.45 


6 


2.65 ± 5.96 


0.26± 0.49 


1.71 ± 1.96 


WEr 


0.575 to 1.96 


0.2 


1.80 


-0.01 ± 0.08 


2.10 ± 0.79 




0.575 to 1.96 


0.9 


1.56± 0.11 


-0.06 ± 0.04 


2.00 ± 0.68 


WE 8 


79 to 2 36 


13 


1.80 


-0.17 ± 0.08 


1.09 ± 0.71 




0.79 to 2.36 


15 


1.52± 5.39 


-0.21± 0.08 


1.08 ± 0.67 


W0 4 


0.89 to 1.46 


23 


1.80 


0.07 ± 0.06 


0.39 ± 0.48 


S* u : WE 5 


0.577 to 1.95 


2.8 


1.80 


0.01 ± 0.01 


0.19 ± 0.13 


WE 8 


0.639 to 1.70 


9.4 


1.80 


0.03 ± 0.01 


-0.18 ± 0.06 


W0 2 


0.846 to 1.38 


18 


1.80 


0.11 ± 0.09 


1.00 ± 0.96 


wo 8 


0.662 to 1.66 


5 


1.80 


-O.Oldb 0.01 


-0.30 ± 0.18 


W0g 


0.627 to 1.73 


3.4 


1.80 


-0.01 ± 0.01 


-0.33 ± 0.19 


E*~: WEi 


0.662 to 1.54 


1 


1.80 


-0.05± 0.19 


-3.34 ± 1.33 


WE 3 


0.926 to 1.42 


16 


1.80 


-0.17± 0.08 


-1.42 ± 0.71 


WE 4 


0.602 to 1.61 


1.3 


1.80 


0.06 ± 0.10 


-1.70 ± 0.38 


WE 5 


0.735 to 1.37 


13 


1.80 


-0.33± 0.36 


-1.40 ± 1.74 


WE 6 


0.588 to 1.97 


0.2 


1.80 


0.01 ± 0.07 


-1.72 ± 0.63 


WE 8 


0.71 to 2.51 


9.4 


1.80 


0.15 ± 0.07 


-1.22 ± 0.65 


WOi 


0.618 to 1.05 


10 


1.80 


-0.34 ± 0.77 


-0.66 ± 1.45 


wo 4 


0.89 to 1.57 


19 


1.80 


-0.08 ± 0.05 


-0.54 ± 0.39 
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TABLE III. 


Same as Table 


| but for 


E*° and E*~. 




bum Kule 


Region 


Cont 


w 


A 






(GeV) 


(%) 


(GeV) 


(GeV x ) 




~* U ; WE 8 


0.636 to 1.55 


13 


2.00 


0.08 ± 0.03 


-0.35 ± 0.12 


wo 2 


0.977 to 1.25 


13 


2.00 


0.11± 0.13 


2.25 ± 1.92 


wo 8 


0.654 to 1.85 


3.7 


2.00 


-0.02± 0.02 


-0.62 ± 0.34 


WOg 


0.621 to 1.91 


2.8 


2.00 


-0.02± 0.02 


-0.69 ± 0.35 


E*~: WEi 


0.628 to 1.61 


1.3 


2.00 


-0.07± 0.14 


-2.88 ± 1.02 


WE 2 


0.898 to 1.12 


37 


2.00 


-0.30± 0.83 


-3.65 ± 5.68 


WE 3 


0.906 to 1.53 


12 


2.00 


-0.12± 0.05 


-1.25 ± 0.55 


WE 4 


0.6 to 1.66 


0.3 


2.00 


0.0004± 0.07 


-1.38 ± 0.27 


WE 5 


0.74 to 1.50 


10 


2.00 


-0.22± 0.21 


-1.27 ± 1.30 


WE 6 


0.59 to 2.11 


0.5 


2.00 


-0.006 ± 0.04 


-1.38 ± 0.48 


WE 8 


0.70 to 2.54 


6.3 


2.00 


0.07 ± 0.05 


-0.88 ± 0.47 


WOi 


0.641 to 1.12 


13 


2.00 


-0.35± 0.57 


-0.58 ± 1.26 


W0 4 


0.863 to 1.86 


12 


2.00 


-0.05± 0.02 


-0.60 ± 0.29 



TABLE IV. Comparisons of decuplet baryon magnetic moments from various calculations: 



this work (QCDSR), lattice QCD (Latt) |15[, chiral perturbation theory (%PT) |L6|, light-cone 
relativistic quark model (RQM) fi~9|| , non-relativistic quark model (NQM) fli~8|| , chiral quark-soliton 



model (xQSM) [24|. All results are in units of nuclear magnetons. 



Baryon 


Exp. 


QCDSR 


Latt 


XPT 


RQM 


NQM 


XQSM 


A++ 


4.5 ± 1.0 


4.13 ± 1.30 


4.91 ± 0.61 


4.0 ± 0.4 


4.76 


5.56 


4.73 


A+ 




2.07 ± 0.65 


2.46 ± 0.31 


2.1 ± 0.2 


2.38 


2.73 


2.19 


A 


« 


0.00 


0.00 


-0.17 ± 0.04 


0.00 


-0.09 


-0.35 


A~ 




-2.07 ± 0.65 


-2.46 ± 0.31 


-2.25 ± 0.25 


-2.38 


-2.92 


-2.90 


£*+ 




2.13 ± 0.82 


2.55 ± 0.26 


2.0 ± 0.2 


1.82 


3.09 


2.52 


s*° 




-0.32 ± 0.15 


0.27 ± 0.05 


-0.07 ± 0.02 


-0.27 


0.27 


-0.08 






-1.66 ± 0.73 


-2.02 ± 0.18 


-2.2 ± 0.2 


-2.36 


-2.56 


-2.69 






-0.69 ± 0.29 


0.46 ± 0.07 


0.1 ± 0.04 


-0.60 


0.63 


0.19 






-1.51± 0.52 


-1.68 ± 0.12 


-2.0 ± 0.2 


-2.41 


-2.2 


-2.48 




-2.024 ± 0.056 


-1.49 ± 0.45 


-1.40 ± 0.10 


input 


-2.48 


-1.84 


-2.27 
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FIG. 1. Diagrams considered for the decuplet baryon magnetic moments. 
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FIG. 2. Diagrams considered for the strange quark mass corrections to the decuplet baryon 
magnetic moments. 
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0.0 0.5 1.0 1.5 2.0 2.5 

M~ 2 (GeV- 2 ) 



FIG. 3. Overlap plots of the valid QCD sum rules for the A ++ magnetic moment. Each sum 
rule is searched independently. The solid line corresponds to the ground state contribution, the 
dotted line the rest of the contributions (OPE minus continuum minus transition). The error 
bars are only shown at the two ends for clarity. From top down, the sum rules are arranged by 
magnitudes of \x B extracted from them. For better viewing, the curves for each sum rule are shifted 
downward by 3 units relatively to the previous one. 
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M~ z (GeV- 2 ) 



FIG. 4. Similar to Fig. §, but for Or. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

M~ z (GeV -2 ) 

FIG. 5. Similar to Fig. |, but for £*+. 
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M~ z (GeV- 2 ) 



FIG. 6. Similar to Fig. g but for 
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FIG. 7. Similar to Fig. |, but for £*~. 
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FIG. 8. Similar to Fig. || but for ~*°. 
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FIG. 9. Similar to Fig. § but for 3*~. 
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FIG. 10. Scatter plots showing correlations between the magnetic moment of and the 
standard QCD input parameters for the sum rule from WE5. The result is drawn from 430 QCD 
parameters sets. 
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FIG. 11. Similar to Fig. [jl| but for S*° and the sum rule from WOs- The result is drawn from 
1000 QCD parameters sets. 
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